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ABSTRACT 

We investigate relativistic flows after a shock wave generated in a star arrives at the surface. First, 
the sphericity effect is involved through a successive approximation procedure by adding correction 
terms to an already known self-similar solution in the ultra-relativistic limit assuming the plane parallel 
geometry. It is found that the lowest order correction terms are of the order of (a/i?*) 1 ^ 6 ^ -8 ) for 
a star with a radiative envelope, where a is the distance from the stellar surface at the moment of 
shock breakout and i?* the stellar radius. We also find that the involvement of sphericity increases 
the acceleration in the early phase as compared with that of the original plane-parallel flow. Second, 
we obtain semi-analytic solutions for a mildly relativistic flow in which the rest mass energy density 
is not negligible in the equation of state. To take into account this, we use the enthalpy and the 
pressure instead of using the density and the pressure as thermodynamic variables. These solutions 
assume self-similar evolutions except for the initial conditions. Third, we have carried out numerical 
calculations with a special-relativistic hydrodynamical code based on the Godunov method in order 
to check the applicability of the above sphericity corrections and the semi-analytic solutions. The 
equation of state used in our calculations includes the rest mass energy density. Comparisons with 
results of numerical calculations support the validity of the sphericity correction terms. The evolutions 
of the pressure and the Lorentz factor of each fluid element of the semi-analytic solution for mildly 
relativistic flows match the numerical results at least in early phases. We also investigate the final 
free expansion phases by this code. For either spherical symmetric or plane-parallel ultra-relativistic 
flows, we could not observe the free expansion phase till the self-similar variables become as large as 
10 6 . On the contrary, the flow in which the ratio of the pressure to density at the shock breakout is 
less than unity reaches the free expansion. We have derived the final energy distributions for these 
flows and compare them with previous works. 

Subject headings: hydrodynamics — shock waves — relativity — self-similar 



1. INTRODUCTION 

Plenty of studies have been made about shock propa- 
gation in stellar explosion. After the energy is liberated 
inside a star, a shock wave is formed to go through to- 
ward the surface. When it reaches there, it disappears. 
This event is called as shock breakout or shock emer- 
gence. Material of the stellar envelope accelerated by 
the shock wave expands into the surrounding space as a 
rarefaction wave. The ejected matters continue to accel- 
erate until their thermal energy is sufficiently converted 
into the kinetic energy. Finally, they expand free ly. 

First. iGandel'man fe Frank-Kamenetskvl (|1956l ) posed 
this problem. ISakuraiT i 1960T ) solved it with a self-similar 
analysis, assuming one dimensional plane-parallel (PP) 
configuration. The stellar surface marks the boundary, 
on one side of it the atmosphere fills and on the other side 
the vacuum. The initial density distribution is assumed 
to follow the power law po oc x , where x is the distance 
from the surface. Especially, 5 — 3 represents the ra- 
diative envelope and S = 1.5 convective. His self-similar 
solution is strictly restricted in applying to a real star, be- 
cause sphericity is no longe r ignorable in describ i ng the 
flow far from the surface. Kazhda n fc Murzinal (|1992t ) 
worked to involve sphericity by adding correction terms 
proportionate to the power of non-dimensional parame- 
ter x = l — r/R*, where r is t he distance from the center 
and R* is the stellar radius. iMatzner fc McKeel (|1999l ) 
conducted a more detailed study. They investigated the 



progenitor dependence and the sphericity involvement 
and derived semi-analytic distributions of physical quan- 
tities such as the velocity and the density. These studies 
ignore effects of the special relativity and the radiative 
cooling, which causes the speed of the front in these stud- 
ie s to be developed to infin ity. 

iTan. Matzner. fc McKeel (|2000t) studied trans- 
relativistic flows. They estimated the maximum velocity 
of ejecta including the effects of sphericity and radia- 
tive diffusion. They also studied effects of aspherical 
structure and gravity. 

In t he ultra-relativistic limit. INakavama fc Shigevamal 

(2005) obtained a self-similar solution ignoring the rest 
mass ener gy. They assum ed the same initial configura- 
tion as in Sakurai (1960). In addition, the equation of 
state (EOS) was assumed as e = 3p in stead of e = 3p+p, 
where p is the pressure, p the density, and e the internal 
energy density per unit volume. This is the in dispensable 
assum ption to derive a self-similar solution. iPan fc Saril 

(2006) presented an analytic expression for the flow after 
shock emergence. 

In this paper, we concentrate on the following two is- 
sues inherent to the above ultra-relativistic self-similar 
solutions. First, the assumption of spherical symmetry 
(SS) is a better approximation of a real system than that 
of PP. The farther a wave is from the stellar surface, the 
larger the deviation from the PP. The end of acceleration 
can not be described within the range that the PP ap- 
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proximation is valid, i.e. |a; | *C 1. The pressure varying 
as p oc p 4 / 3 oc (1 — iTo) -4 does not significantly drop in 
that range. To take into account the sphericity, we need 
to introduce a characteristic length, the stellar radius i?*, 
which inhibits analyses by self-similar solutions. Thus, 
in <J3J we conduct a successive approximation in terms of 
a parameter £ = a/i?* ( where a is the distance of each 
fluid element from the stellar surface at shock emergence 
) by which the sphericity is characterized. From the set 
of fundamental equations for ultra-relativistic hydrody- 
namics, we found that the sphericity is involved from the 
order 0(^/((2V3-3)«+i)) after ghock ernergence anc ; 

before shock emergence. We take this effect into account 
only up to the order 0(^/((2V3-3)«+i))_ 

Second, the acceleration does not terminate in 
iNakavama &; Shigevamal (|2005f) 's solution. The Lorentz 
factor 7 of a fluid element increases infinitely with in- 
creasing time t measured from the moment of shock 
emergence; 7 oc i - 732 for any S. The source energy of 
acceleration is the internal energy which is converted to 
the kinetic energy during expansion. It is reasonable to 
consider that when the internal energy is sufficiently con- 
sumed, the ejecta are no longer accelerated. This implies 
that their self-similar solutions can not describe a stage 
of free expansion. It might result from the EOS e = 3p 
neglecting the rest mass energy. The EOS has to be 
£ = 3p + p in later phases when the gas becomes mildly 
relativistic. In Sj3l we derive the semi-analytic solution 
including this term in the EOS . 

In <21 we conduct numerical calculations using the Go- 
dunov method for relativistic hydrodynamics. The aim is 
to support the validity of the sphericity involvement and 
the semi-analytic solutions for mildly relativistic flows. 
Calculations are carried out from the moment of shock 
emergence for both PP and SS configuration and vari- 
ous initial energy strengths, setting S = 3 and the ini- 
tial conditions from the analyt ic solutions at t = of 
INakavama fc Shigevamal ((2005) . The results are pre- 
sented in 



2. SPHERICITY INVOLVEMENT 
2.1. Formulation 

In this section, we derive sphericity corrections to the 
PP flow after shock emergence. First, we present the 
basic equations for an ultra-relativistic SS flow. Sec- 
ond, we transform these equations into non-dimensional 
ones. Third, we conduct a successive approximation 
to obtain 0th and 1st order equations. The 0th order 
equations provide the se l f-simi lar solutions obtained by 
INakavama fc Shigevamal (|2005f K The 1st ones shall pro- 
vide the amount of corrections. 

We confine ourselves to dealing with a perfect fluid 
with no external force including gravity. If a shock wave 
is strong enough, gravity has practically no effect on the 
fluid motion. We assume the initial density distribution 
po = b(R*—r) s . The basic equations of relativistic hydro- 
dynamics are derived from the energy-momentum and 
mass conservation. Here we present the equations for 
spherically symmetric flows in a fixed frame (r, t) . 



0_ 

di 



1 d 



[ 1 2(e + f3 2 p)} + -—[r 2 1 2 (e+p)p}=0, (1) 



I h 2 (e +P )p] [rVCe + ri/J 8 ] + % = 0, (2) 



~dt 



(3) 



where e denotes the internal energy density, p the pres- 
sure and p the density, all measured in the fluid frame. 
p' = 7p is the density measured in the fixed frame where 
the fluid moves with the velocity c/3. 7 denotes the 
Lorentz factor defined as 7 = - (3 2 . We set the 

speed of light c to 1 in the following. We assume the 
ultra-relativistic EOS as 



3p. 



(4) 



With this EOS (gj), equations (dJ, and © 

are transformed in t o th e expressions presented in 
iBlandford fc McKed (fl976h : 



dt^ 



2 dp 

7 at' 



4MpV) = -~|V/3), 



dt 



dr 



dt 



(pp" 4/3 ) = o, 



where the convective derivative 

— - — +13— 
dt dt dr ' 



(5) 



(6) 



(7) 



(8) 



has been introduced. 

In order to describe the motion after shock emergence, 
it is convenient to use a Lagrangian coordinate a which is 
the position of a fluid element at t = measured from the 
stellar surface. Then from the definition of the velocity 
of the fluid clement, we obtain 



dr 
di 



(9) 



For a PP flow, only one non-dimensional variable rj is 
sufficient to describe the self-similar motion of the gas. 
In contrast, for an SS flow, when we deal with the wave 
across the stellar surface, the stellar radius i?* appears 
as a scale. The flow is not self-similar anymore. Another 
non-dimensional variable is necessary besides ry. We can 
write the non-dimensional variables such as 

v=4-> ( io ) 



7fa 
a 



(11) 



The dependent variables are written using non- 
dimensional functions, F, G, H and /, such as 



p = Pi(a)F(ri,g), 

7 2 = 7i 2 («)£(??, a 

p = pi (a)H (r),£), 
r=(R*-a)I{ri,0. 



(12) 
(13) 
(14) 
(15) 



Here Pi(a), 7 2 (a) and pi(a) are the profiles at t = 0. 
They conform to the power law, 



p l {a)=c p {A^W^ m)/{m+1 \ 



(16) 
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1 2 1 (a) = c 1 (A,b)a- m ^ m+1 \ 



(17) 



p[ (a) = 71P1 = c p , (A, b ) a ^-m)/{m+i) , (18) 

where 

m= (2V3-3)S, (19) 

and c p (A, b), c 7 (A, b) and c p ' (A, b) are constants depend- 
ing on A and 6. These two quantities characterize the 
flow before shock emergence. The parameter A speci- 
fies the amount of energy of the flow, through a relation 
with the shock Lorentz factor f by f 2 = A(—t) m . For 
more detailed descriptions, see Nak avama fc Shigevamal 
(2005). We transform the independent variables from 
(i, a) to (?7,f) using the following relations, 



d d drj 3 <9£ ry d 
It = d~q~di + ~dl~di = J drj' 



d d drj d <9£ 



1 r) 8 £ d 



da drj da d^ da m + 1 a drj a <9£ 



(20) 



(21) 



Then equations ([5j) , ([6]) , (JTJ) and ((9]) are converted to the 
following non-dimensional equations. 



— (lnF + 21nG) = 
drj 



I 2 H /8-m rj dF dF 
\m + l m+ldr/ d£ 



, (22) 



<91nG 3<91nF 



drj 2 drj 

\ m + 1 drj dt; m + 1 

4 7l 2 e g 



<91nF AdlnH 



drj 3 drj 

01 _ lit 
drj 1-f 



, (23) 



(24) 



(25) 



The expansions in 7~ 2 can be truncated at the first 
contributing term because the motion is assumed to be 
ultra-relativistic. The last term of (|23|) is of the order 
Q{^/ ("*+!)). Therefore we expand the solutions (Jl2j) - 
(|15|) in a power series of £ starting with the term propor- 
tional to 

F( V ,0=F a ( V )(l+e /(m+1) F 1 ( V )), (26) 

G(r/,a-G (ry)(l+e 1/(m+1) G 1 (r,)), (27) 

H(fl,Z) = H ( V )(1 +e /(m+1) Hx(v)), (28) 

7(r,,e)-/oW(l+e 1/(m+1) /iW). (29) 
The 0th order equations become 
dl 



drj 



0, 



(30) 



dlnF , din Go 



(/;/ 



V dF 
Go V m + 1 m + 1 dry 



(31) 



din G 3dlnF 



dry 



2 dry 



o-noV^o / dlnGo 

' rj — ; m 



m + 1 ^ ' dry 
dlnF 4dlniJ 



(32) 



(33) 



dry 3 dry 
These correspond to the planar ones. The 1st order are 



dh _ m /{m+l) 



(34) 



(m+l)v^ 



dFi 
dry 



,dGi 
dry 



(2/ 1 +lf 1 --G 1 )U-m-ry 



dFo 
dry 



dF dF : 

^ vFo-r- 

drj drj 



FqFi, (35) 



'Gq 



dGi 
dry 



3 dFj 
2~<kj~ + 
dGi 



4c 7 i?; 



-m/ (m+l) 



Go 



dry 



Gi 



1 



(2h +H 1 + -G 1 



dlnGo 
dry 



dFi 
drj 



AdH x 
3 dry 



(36) 



(37) 



The solutions of (j30|) - (f37|) are obtained through numeri- 
cal integration with the following initial conditions; 

F (0) = G o (0) - ff (0) = /o(0) = 1, (38) 
Fi (0) = Gi (0) = Hx (0) = 7i(0) = 0. (39) 

This implies that we adopt the distributions of PP ge- 
ometry at t = as the initial condition of SS geome- 
try, because we can neglect deviations from PP flows at 
shock emergence. All sphericity corrections before shock 
emergence are 0(£) that is higher than that after shock 
emergence, 0(£V("H-i)). 

2.2. Solutions 

The results are shown in Figure [1] The dash-dotted 
lines represent the flow involving the sphericity correc- 
tion, and the dashed lines the self-similar solutions. The 
density and the pressure involving the effect of sphericity 
drop more steeply than in PP. 

It is appropriate to confirm the range in which the 
correction is valid. This is one of incentives to conduct 
numerical calculations. In ij4.2.2[ the validity of our cor- 
rection will be examined. 
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Fig. 1. — Comparison between a planar flow and a spheri- 
cal flow. The dashed, the dash-dotted, the solid and the dotted 
lines correspond to the self-similar solutions, the solutions involv- 
ing sphericity, an ultra-relativistic PP flow from a numerical cal- 
culation and, an ultra-relativistic SS flow from a numerical calcu- 
lation, respectively. These two numerical calculations were carried 
out with the parameter A = 1.0 X 10 16 . The solid and dotted 
lines from the results of the numerical calculations for PP or SS 
are drawn by tracing one fluid element with the preshock position 
xo/Rt = 4.49 X 10 -4 and at shock emergence pi/pi = 33.1. The 
sphericity correction is obtained for £ = 1.15 X 10 — 6 . 



3. MILDLY RELATIVISTIC FLOW 

3.1. Formulation 

In later phases, p becomes comparable to p. We will 
not be able to ignore the rest mass energy density in 
the EOS. We should use the EOS e — 3p + p instead of 
s = 3p. In order to include the term of the rest mass 
energy density, we replace the dependent variable p in 
equations (P)-© with the enthalpy per unit volume : 
w = Ap+p. As a result, the basic equations are converted 
to 

: (w 1 2 (3) + —(w 1 2 (3 2 +p) = 0, 



d_ 

di 



d 



{wY - p) + —{wY P) =o 



dx 



d d 
— [7O - 4p)] + 



Ap)} = 0. 



These equations are combined to yield 

1 dp _ dp 

dx ' 



w_drf_ 2 dp 
2 dt 7 dt 



1 dhi7 2 
2 



dlrtw 



dt 



dt 



w dt 



4*>)] + -4(p- 
w dt 



w) 



0. 



(40) 
(41) 
(42) 

(43) 
(44) 
(45) 



To use the Lagrangian coordinate as in the self-similar 
ultra-relativistic PP solution, another equation not con- 
taining p is needed. This equation can be obtained by 
differentiating the definition of the fluid velocity with re- 
spect a as 

d dx dv , 

(46) 



dt da da 
We seek a solution in the form 

P = Px(a)F{rj), 

7 2 =7i 2 («)G(r/), 

w = 4pi(a)W(r)) 

dx 

l)a 

by analogy with the ultra-relativistic case. We substi 



(47) 
(48) 
(49) 

(50) 



tute 1|17 |) -([5D ]) into 
equations, 

2W G' 
~F~~G 



and obtain non-dimensional 



1 



(m + l)GJ 



2 1 



V 



1 



IGJJ G 
W 



W 
W 



EL 

F 

EL - 
w ~ 

EL] EL 

WJ ~F 



1 



m+lGJ' 
2m 1 



IGJ 



J' = 



1 



2(m+ 1)G 
The initial conditions are 

F(0) = 1, 
G(0) = 1, 

W(0) = 1 

J(0) = 1. 



G' 
V G 



4pi(a)' 



(51) 
(52) 
(53) 
(54) 

(55) 
(56) 

(57) 

(58) 
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TABLE 1 

Power-law indices of F, G, W and J for the semi-analytic 

SOLUTIONS 



10" 1 10° 10 1 10 2 10 3 10 4 10 5 10° 



FlG. 2. — Comparison between the semi-analytic solutions for 
eqs. H5U - H54H and the self-similar solutions. The solid, the dotted, 
the dash-dotted, and the dashed lines represent the semi-analytic 
solutions with the initial conditions pi/ pi = oo, 2.5, 0.5, 0.25 re- 
spectively. The crosses, +, represent the self-similar solutions. 



If W(0) = 1, the solutions agree with 

iNakavama fc Shigevamal {2005). Otherwise, they 
are not self-similar and F, G, W, and J depend not only 
on rj but also on a and A through the initial condi- 
tion (|57p . We will check the validity of the assumed 
form of the solution by comparison with results ob- 
tained from direct numerical integrations of rclativistic 
hydrodynamical equations in 

3.2. Solutions 

F and G integrated with various W(0) are shown in 
Figure [2] It is obvious that mildly relativistic flows, un- 
like an ultra-relativistic flow that forgets its history, de- 
pend on the initial energy distributions. 

From the numerical integration, it is found that 
F, G, W and J undergo power-law evolution in the limit 
of rj — > oo. The indices are shown in Table [1] for rj > 10 6 . 
The values of d In G/d In r/ in this table indicate that even 
though the semi-analytic solutions include the rest mass 
energy in the EOS, the acceleration does not terminate 
completely. For pi/pi 1, the acceleration almost ter- 



Pi/ pi 


d\nF/dlnri 


din G/d ln V 


d\nW/d\nr) 


d In J/d In rj 


oo 


-0.845 


0.732 


-0.845 


0.268 


2.5 


-1.01 


0.304 


-0.813 


0.607 


1.0 


-1.13 


0.159 


-0.8744 


0.772 


0.5 


-1.24 


0.07 


-0.935 


0.893 


0.25 


-1.27 


0.04 


-0.957 


0.931 



mutates for rj — > oo. In contrast, the acceleration does 



not terminate for pi / pi 



1 



dV 


dv 


dt ~ 


dm x ' 


ds 


dp 


~dl~ 


dm x ' 


dQ _ 


d{pv) 


dt ~~ 


dm x 



4. NUMERICAL CALCULATIONS 

4.1. Calculation procedure 

We carry out numerical calculations with the Go- 
dunov method for relativisti c hydrodynamics in th e La- 
grangian coordinate system (|Marti fc Mu llcr 1995j). To 
write down the basic equations for relativistic hydrody- 
namics in conservation form, we have to introduce some 
new variables V, s, and Q. These variables are related 
to quantities p and w in the local rest frame of the fluid 
through 

V=l/(p 7 ), (59) 
s = (w~fv)/p, (60) 

Q = ( W1 2 -p)/(p7)-l. (61) 

For PP geometry, the basic equations are expressed with 
these variables as 

(62) 
(63) 
(64) 

(65) 

(66) 
(67) 
(68) 

(69) 

These equations are integrated numerically by evaluating 
the derivatives in the right hand side using the Godunov 
method. The difference lies in the EOS; e = 3p + p in 
numerical calculations while e = 3p in the self-similar 
solutions. 

Calculations were carried out for both PP and SS ge- 
ometry, from the moment of shock emergence (t = 0). 
We set the initial condition from the self-similar solutions 
at the time t — 0. The parameters used for calculations 
are 5 = 3, A = 1.0 x 10 16 , 1.0 x 10 12 or 1.0 x 10 10 , and 



where the Lagrangian coordinate is defined as 

dm x = — pi^/ida = —pjdx. 
For SS geometry, 



dV 




~dt~ 


/ <* 


ds 


- 9 


~dl~ 


dm r 


dQ _ 


- 9 


dt ~ 


dm r 



(47rr 2 p), 
(47rr 2 pi>), 

where the Lagrangian coordinate is defined as 
dm r = —4it(R* — a) 2 pi"fida = Airr 2 p-fdr. 
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i?* = 1.0 x 10 9 m. Here A has a dimension of some power 
of time. Thus the values of A are expressed using second 
as a unit of time. Note that 71 depends on A as 



7i oc 



A" 



(70) 



Hence 71 varies by a factor of about 10 for the above 
values of A. 

We obtain p, p and 7 directly from numerical results. 
The non-dimensional variables F, G, and H can be ob- 
tained by p/pi, 7/71, and p/pi, respectively. On the 
other hand, when we made calculations starting before 
shock emergence, the moment of shock emergence could 
not be determined exactly. As a result we could not de- 
termine pi, 71, or pi to sufficient accuracy. This was 
caused by the finite thickness of the shock front in nu- 
merical calculations. Therefore, we dedicate ourselves to 
calculations after shock emergence. 

4.2. Results 
4.2 .1. Connection to the initial states 

In presenting the results after shock emergence, it is 
preferable to describe those in terms of the original posi- 
tion of each fluid element xq rather than the Lagrangian 
coordinate a to help actual applications to astrophysi- 
cal phenomena. For PP, the definition of the Lagrangian 
coordinate provides the exact relation between a and xq. 



padxa 



p'ida, 



(71) 



Substituting po = bx s and the expression for p' x (Eq. 
(|18p) into both sides of this equation, we obtain 

x Q = [2A(m+l)] 1/{S+1) 

( 2A[1 + (2^3 - 3)6] 



(m+l)(« + l) 



xa 



[3 + 2^3(1 + ( S)]v / 3(l+<5)/{2[l+v / 3(l+«5)]} 
l/(m+l) 



If S = 3, then x = 1.56 x A 



418 



(72) 

x a 418 . Though the 



corresponding relation ro = ro (a) for SS will be obtained 
from 



I 4tt(R« - a) 2 p[da = / 471Tq 

J0 Jra 



p dr , (73) 



it is difficult to express ro as a function of a in an analytic 
form unless we truncate higher terms of a/ R„ and Xq/ R*. 
To the lowest order, it takes the same form as PP, 



r Q ~[2A{m+l)} 1/{6+1) 

( 2A[l+ (2V3-3)<5] 



(m+l)(S+l) 



xa 



[3 + 2\/3(l + ( 5)]v / 3(l+<5)/{2[l+V^(l+<5)]} 
l/(m+l) 



(74) 



Using these relations, a fluid element after shock emer- 
gence can be traced back to its original position. 

In addition, we can see how the parameters character- 
izing the explosion such as the injected energy E- ln j and 
the ejected mass M e j C affect the flow through the param- 
eter A as follows. The parameter A can be related to 



TABLE 2 

Differences between the numerical 
solutions for a spherical flow and the 
self-similar solutions with correction 
terms for sphericity 



V 


AF/F (%) 


AG/G (%) 


AH/H (%) 


0.3 


2.1 


1.0 


1.2 


0.4 


3.9 


1.8 


2.4 


0.5 


6.1 


2.5 


3.9 



.Einj and M ejc by 



mj = / 7 pdx, 

ocr 2 (-t)-"' 

oc AM eic . 



CJC, 



Thus we obtain 



.4 



M fi 



(75) 



(76) 



It is desirable to get the proportional constant of (|T6"j) 
which connects the shock velo city with E- m ^ and M e j c . 
iTan. Matzner. k. McKeei ([2000) proposed an analytical 
formula (eqns (4) and (9) in their paper) to estimate 
such a coefficient. See a ls o equ ations (17), (18) and 
(19) in iMatzner fc McKed (|1999h . In our case, we can 
obtain an expression of A as a function of E^ and 
M e j c by substituting p = b(— t) s into equation (4) in 
ITan. Matzner. fc McKeei pOOO ). 

4.2.2. Effects of sphericity 

We discuss the validity of our sphericity correction in 
this section. Figure [1] shows the trajectories of the fluid 
elements with the preshock position xq/R* = 4.49 x 10 -4 
and at shock emergence pi/ pi = 33.1 (The ratio of the 
internal energy to the kinetic energy is about 132). The 
solutions of (|2"u]) - ([2"8"j) including sphericity correction (the 
dash-dotted line) match the numerical results for SS with 
A = 1.0 x 10 16 (the dotted line) for small 77. On the con- 
trary, the planar flows start to deviate from the begin- 
ning. The degree of accuracy in the sphericity correction 
is shown in Table [2l Our description for an SS flow is a 
good approximation in the range of 77 £ 0.5 within accu- 
racy of 5 % . 

Next, we compare an ultra-relativistic PP and SS flow 
from the results of numerical calculations shown in Fig- 
ured] The density of SS decreases more drastically than 
that in PP (the solid lines). This results from the dif- 
ference in the geometry; the volume of a spherical shell 
expands more than that in a PP shell. Then, the pres- 
sure in SS also decreases more than that in PP, because 
the flow is adiabatic. The Lorentz factor in SS increases 
more enormously than that in PP. As the acceleration is 
advanced, the ratio of the internal energy to the kinetic 
energy which is about 4p/p for 7 — * 00 in the observ- 
ing frame drops with decreasing p because pj p decreases 
as p 1 / 3 in an adiabatic flow. As a consequence, the ac- 
celeration of an SS flow is more quickly advanced and 
terminates in a shorter time scale than PP. 

4.2.3. Intermediate distribution 
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Fig. 3. — Numerical results of an ultra-relativistic PP flow for 
A = 1.0 X 10 16 . Time evolutions of F, G and H are shown. Each 
line represents the locus of one fluid elements with the preshock 
position from xq/R* = 6.27 X 10~ 4 to 1.36 X 10~ 2 and at shock 
emergence pi /pi ranging from 3.09 to 26.3. A considerable overlap 
of all the lines suggests self-similarity of the flow. 

For r] — ► oo, F, G, and H from the self-similar so- 
lution have power law distributions with respect to r). 
We call such distributions as the intermediate distribu- 
tions. F, G, and H from numerical calculations can be 
fitted by power law distributions. The indices are shown 
in Table [3] which are obtained for 77 ^ 10 — 100. Com- 
paring the self-similar flow and the ultra-relativistic PP 
flow obtained from numerical calculations, the differences 
of indices between them are 1.1%, 5.7%, and 0.63% for 
dlnF/dlnr/, dlnG/dlnr], and d In H/d In 77, respectively. 
The difference in d In G/d In rj is particularly prominent, 
which implies that the effect of the rest mass energy den- 
sity shows up to suppress the acceleration. That is to 
say, d In G/d In 77 of numerical calculations is smaller than 
that of the self-similar solutions. The values of Pi/pi 
ranges form 10 to 70 depending on a in this particular 
numerical calculation. 

Figure[3]shows the evolution of several fluid elements in 
a PP flow obtained from a numerical calculation. Though 
each line has a different value of pi/ pi in the range from 
8.5 to 33, the trajectories are insensitive to it, which in- 
dicates the flow is very close to self-similar. If we set 
Pi/ Pi ~> °°: the result will approach the self-similar so- 
lutions. 

Numerical results of an SS flow behave differently from 
that of a PP flow. These are shown in Figure HI The 
trajectory of each fluid element of an SS flow depends 
on the initial position a. It is insufficient to use 77 for 
characterizing the flow. The values of indices change in 
the ranges shown in Table [3] 

4.2.4. Mildly relativistic flow 

The acceleration in numerical calculations is sup- 
pressed compared to that of the self-similar solutions due 
to the rest mass energy in the EOS used in numerical cal- 




FlG. 4. — Same as Figure [3] but for a UR spherical flow with the 
preshock position from xq/R* = 4.49 X 10 -4 to 3.29 X 10~ 3 and 
at shock emergence pi/pi ranging from 8.48 to 33.1. 



TABLE 3 

Power-law indices of F, G, and H for ultra-relativistic 

FLOWS IN THE INTERMEDIATE ASYMPTOTIC PHASE 





din Fj dint] 


din G/d In 77 


d\nH/d\nr) 


self-similar 


-0.845 


0.732 


-0.634 


PP (numerical) 


-0.836 


0.690 


-0.630 


SS (numerical) 


-3.23 - -2.97 


1.33 - 1.44 


-2.43 - -2.27 
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Fig. 5. — Numerical results for a mildly relativistic PP flow. 

This calculation was carried out with A = 1 X 10 10 . The solid, the Fig. 6.— Numerical results for a mildly relativistic SS flow, 

dotted, and the dash-dotted lines correspond to the fluid element This calculation was carried out with A = 1 X 10 1(1 . The solid, the 

with the preshock initial position, x /R* = 1.05 X 10 3 , 1.72 X dotted, and the dash-dotted lines correspond to a fluid element 

10 -3 , 2.88 X 10 -3 , and at shock emergence the strength of rela- with the preshock initial position, xo/R* = 7.96 X 10~ 4 , 1.33 X 

tivity, pi/pi = 1.022, 0.727, 0.507, respectively. i0 - 3 , 2 .37 x 10" 3 , and at shock emergence the strength of rela- 
tivity, pi/pi = 1.24, 0.870, 0.580, respectively. 



TABLE 4 

Power-law indices of F, G, and H for mildly relativistic flows 

IN THE FINAL FREE EXPANSION PHASE (FROM NUMERICAL RESULTS) 





xq/R* 


Pi/ Pi 


dlnF/dlnn 


dlnG/dlnn 


din H/dlnn 


PP 


1.05 x 10" 3 


1.02 


-1.21 


0.03 


-0.98 




1.72 x 10" 3 


0.727 


-1.25 


0.02 


-0.98 




2.88 x 10" 3 


0.507 


-1.26 


0.03 


-0.97 


SS 


7.96 x 10" 4 


1.24 


-3.31 


0.00 


-2.99 




1.54 x 10~ 3 


0.782 


-3.45 


0.00 


-2.99 




2.59 x 10~ 3 


0.546 


-3.52 


0.00 


-3.00 



4.2.5. Final distribution 



We investigate the final energy distributions for the 
mildly relativistic flow and the PP flow from numerical 
calculations. The energy per unit mass excluding the 
rest-mass energy is given by 



where 



PI 

t = (p + 4p)7 2 — p - fry. 



(77) 
(78) 



culations. 

The results of mildly relativistic flows from numerical 
calculations are shown in Figure [5] for a PP configura- 
tion and Figure [5] for an SS configuration. In Figures 
[7] and we compare the semi-analytic solutions (PP) 
for a mildly relativistic flow with the corresponding re- 
sults from numerical calculations. The power law indices 
in the intermediate distribution are shown in Table [4) 
The semi-analytic solutions give a precise description of 
a mildly relativistic flow up to rj ~ 100. 



In the stage of free expansion, the ejecta are accelerated 
to 7 = 7f >■ 1. Thus, 7f ~ e. Using the definition of the 
Lagrangian coordinate, we calculate the energy spectrum 
defined as the total mass M(> 7f) with the energy per 
unit mass greater than 7f ~ e as 



M(> 7f ) = / pjtdx, 



2.6 



(79) 



Furthermore, it is useful to present the energy of the 
ejecta E(> 7) with Lorentz factors greater than 7. From 
the numerical results for a PP mildly relativistic flow, 
when 77 — > 00 we obtain the relation between a and the 
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Fig. 7. — Comparison between the self-similar solutions, the 
semi-analytic solutions for a mildly relativistic PP flow, and the 
numerical results for a mildly relativistic PP flow. The dotted 
lines, the solid lines and the crosses, +, represent the self-similar 
solutions, the semi-analytic solutions, and the numerical results, 
respectively. This numerical calculation was carried out with A = 
l.Ox 10 10 . And the crosses represents atrajectory of a fluid element 
with x /R* = 0.00325 and pi/pi = 0.466. 



Fig. 8. — Same as Figure[7]but for a more extended range of r). 



(80) 



final Lorentz factor 7f, 

7f oc a" ' 64 . 
Using this relation, we obtain 

E (> 7) = / UPlfdx, 



2+(m + l)5 

o^hA 7" 1 - 6 . (81) 

iTan. Matzner. fe McKed (|2000f ) derived an empirical 
formula for the final Lorentz factor 7f of each element as a 
function of important parameters such as poj ^inj , -^eje, x 
that characterize the shock emergence phenomenon. For 
a PP flow, They obtained 

7f ~ V 1 + it 2 - 03 + M 1 + 9 2 )°' 12 ) VI ]'?(i + g 2 ) - 12 } 2 , 



where q is defined as 



q = 0.736 




mM e _ 
Po(R*-x) 3 



-1 0.187 



(82) 



(83) 



Here m(r) is the fraction of M e j c 
remnant mass M rcm , i.e. rh(r) 



within r excluding any 

= [m(r) - M rem ]/M e jc 



£"i n j is the injected energy in units of the rest 



mass energy of M cjc , i.e. £ inj = £^ inj /M cjc . In 
the limit of q — > 00, thi s equation tends to reproduce 
Uohnson fc McKed (|197lD and thus M{> 7f) cx jf xl , 
which is also derived from an ultra-relativ istic self-similar 
solution (Nakayama & Shigeyama 2005). On the other 
hand, Sakurai (1960)'s solutions give M(> 7f) oc 7 f ~ 3 ' 6 . 
These results suggest that more relativistic the shock is, 
more energies can be transferred to the matter closer to 
the surf ace. This fact is already pointed out and formu- 
lated bv lTan. Matzner. fc McKeei (l2000h. See their equa- 
tions , (35), (36), and (37), (see also Uohnson fc McKed 

Eml). 

Next we discuss the energy spectrum for spherical 
flows. From the numerical result (A = 1 x 10 10 ) for 
the mildly relativistic SS flow, 



7f tx a 



-0.62 



(84) 
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Fig. 9. — The deficit factors / ap h (see text for detail). These 
results are obtained by numerical calculations for both PP and SS 
configurations with A = 1.0 X 10 10 (solid line) and A = 1.0 X 10 12 
(dotted line). 



The energy spectrum for SS is 

M(>7/)= J p"ff4nr 2 dr, 

= J /3i7i47r(i?* — a) 2 da, 



2.7 



(85) 



~47ri? 2 / pijida, 

ocM (^+i)/( m +i) 7 - 

In terms of the energy, 

E(> 7 ) ex j- 1 - 7 . (86) 

Although the definition of M(> jf) in the PP con- 
figuration differs from that in the SS configuration, the 
difference in geometry is negligible as far as we deal with 
very thin shells. Sphericity does not significantly alter 
the exponents of 7f in E(> 7f) and M{> 7f ) from those 
of the PP flow, but reduces the absolute values of E(> ) 
and M(> 7f ) compared to the PP flo w. 

Following (Tan. Matzner, fc McKee] (|2000l) . we define 
/ sp h which is a deficit factor to characterize the reduction 
of the final velocity of an SS flow from the corresponding 
PP flow, 

spherical (Q'7\ 

Jsph=—, ■ l°'J 

(7fPf Jplanar 

We show the deficit factor / sp h of the flows with A = 
1.0 x 10 10 and A = 1.0 x 10 12 in Figure H This figure 



clearly shows that the effect of sphericity is more promi- 
nent in the outer layers, which was seen in the flows with 
E inj = 0.3 calculated bv lTan. Matzner. fc McKei (|2000h . 
On the contrary, they also showed that non-relativistic 
SS flows become identical to the corresponding PP flows 
in the outer layers. This figure also indicates that the 
spherical expansion in a more intense explosions with 
larger A suppress the acceleration to a greater degree. 
Note, however, that the PP flow with A = 1.0 x 10 12 
does not reach the free expansion phase in our calcula- 
tion. Thus the deficit factor f sp ^ for A = 1.0 x 10 12 shown 
in this figure will be still decreasing as time elapses. 



5. CONCLUSIONS AND DISCUSSION 

We have derived the formula for sphericity correction 
to ultra-relativistic PP flows and presented the semi- 
analytic solutions for mildly relativistic flows after shock 
emergence. Both are proved to have a good agreement 
with numerical results at least in the earliest phases of 
evolution after shock emergence. 

Although the semi-analytic solutions include the rest 
mass energy in the EOS, this inclusion of the rest mass 
energy does not lead to terminating the acceleration if 
the pressure to rest mass energy density ratio pi/pi is 
large enough. The semi-analytic solutions can not de- 
scribe the whole evolution of a flow for finite values of 
Pi/pi. This may be attributed to the assumption that 
a flow is described with one non-dimensional parame- 
ter except the initial condition. Reaching the stage of 
free expansion, a PP flow appears to deviate from self- 
similarity even if it is ultra-relativistic. Furthermore we 
neglect the effect of radiative diffusion which plays the 
critical role to suppress the acceleration. 

We conducted the numerical calculation for an ultra- 
relativistic flow by setting A = 1.0 x 10 12 . The calcula- 
tion continued till r\ reached ~ 10 6 , but it could not reach 
the stage of free expansion. We could not confirm by 
numerical calculations whether an ultra-relativistic flow 
leads to free expansion or not. 

On the other hand, when a flow is not so ultra- 
relativistic at the time of shock emergence, the accel- 
eration comes to an end both in SS and in PP configura- 
tion. However there remains a difference between them; 
a spherical flow rather suddenly finishes the accelera- 
tion while a PP flow approaches the end of acceleration 
asymptotically. 
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